ROSAT HRI Detection of the 16 ms Pulsar PSR J0537-6910 Inside SNR N157B 
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ABSTRACT 

Based on a deep ROSAT HRI observation, we have detected a pulsed signal in the 
0.1-2 keV band from PSR J0537-6910 — the recently discovered pulsar associated with 
the supernova remnant N157B in the Large Magellanic Cloud. The measured pulse 
period 0.01611548182 ms ± 0.02(ns), Epoch MJD 50540.5, gives a revised linear spin- 
down rate of 5.1271 x 10~ 14 s s _1 , slightly greater than the previously derived value. The 
narrow pulse shape (FWHM ~ 10% duty cycle) in the ROSAT band resembles those 
seen in both RXTE and ASCA data 2 keV), but there is also marginal evidence 
for an interpulse. This ROSAT detection enables us to locate the pulsar at R.A., Dec 
(J2000) = 5' l 37 m 47 s .2, -69 o 10'23". With its uncertainty ~ 3", this position coincides 
with the centroid of a compact X-ray source. But the pulsed emission accounts for 
only ~ 10% of the source luminosity ~ 2 x 10 36 ergs -1 in the 0.1-2 keV band. These 
results support our previous suggestions: (1) The pulsar is moving at a high velocity 
(~ 10 3 km s" 1 ); (2) A bow shock, formed around the pulsar, is responsible for most 
of the X-ray emission from the source; (3) A collimated outflow from the bow shock 
region powers a pulsar wind nebula that accounts for an elongated non-thermal radio 
and X-ray feature to the northwest of the pulsar. 



Subject headings: pulsars: general — pulsars: individual (PSR J0537-6910) — X-rays: 
general — supernova remnant 



1 



1. Introduction 

PSR J0537-6910 is the most rapidly rotating 
young pulsar known. It was detected serendipi- 
tously by Marshall et al. (1998) while searching 
for a pulsed signal from nearby SN 1987A with 
the Rossi X-ray Timing Explorer (RXTE). Using 
archival data of the region acquired with the Ad- 
vanced Satellite for Cosmology and Astrophysics 
(ASCA), they further located the pulsed emis- 
sion to within the ~ 1' radius of the supernova 
remnant (SNR) N157B in the Large Magellanic 
Cloud (LMC). This SNR is a moderately bright 
X-ray source in the 30 Doradus star formation 
region and is clearly imaged with the ROSAT 
High Resolution Imager (RHRI) as isolated dif- 
fuse emission (Fig. 1), ~ 14' east of SN 1987A, 
and ~ 15' northwest of the PSR B0540-69. 

Before the discovery of PSR J0537-6910, Wang 
k Gotthelf (1998) showed that the SNR N157B 
contains a bright, elongated, and non-thermal X- 
ray feature, which dominates the SNR emission 
and is most likely a pulsar wind nebula. Near 
the southeast end of the feature is a compact X- 
ray source with a spatial extent of <^ 7". They 
suggested that this source represents the pulsar 
and its interaction with the surrounding medium. 
The predicted SNR age 5x 10 3 yrs agrees with the 
characteristic spin-down age of PSR J0537-6910 
(Marshall et al. 1998). These studies have es- 
tablished N157B as the latest addition to the 
rare class of Crab-like remnants (e.g., Seward 
1989), of which only three confirmed examples 
were known previously: the Crab Nebula, the 
MSH 15-52 nebula, and LMC SNR B0540-693. 

We report in this Letter the detection of pulsed 
soft X-ray emission from PSR J0537-6910, using 
a deep RHRI observation recently acquired. This 
detection is the first in the energy band below 
2 keV. This band-pass is not available to X-ray 
detectors with poor spatial resolution, because of 
the overwhelming soft X-ray emission from both 
the pulsar wind nebula and the diffuse thermal 
hot gas of N157B. Pulse phase-resolved high res- 
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Fig. 1. — Coadded RHRI image of N157B. The image 
is adaptively smoothed with a Gaussian of adjustable size 
to achieve a constant local count-to-noise ratio of ~ 6 over 
the field. Each contour is 2a above its previous level. The 
lowest level is about 2a above the local background of 
~ 7 x 1CP 3 counts s _1 arcmin~ 2 , while the highest one is 
at 1.06 counts s~ x arcmin~ 2 . 

olution imaging capability of the RHRI further 
enables us to locate the pulsar with a linear po- 
sition accuracy of <J 3", a factor of > 10 better 
than the previous measurements. The pulsar lo- 
cation confirms our early suggestion of the com- 
pact X-ray source as the site of the pulsed emis- 
sion. The accuracy of the location may also allow 
for identification of pulsar counterparts in other 
wavelength bands. 

2. Observation and Data Reduction 

The new RHRI observation was taken as part 
of our ongoing X-ray investigation of the 30 
Dor region. The high resolution X-ray image 
of N157B as shown in Fig. 1 is a co-add of 
the data from ROSAT observation Sequence No. 
rh600228n00 (30 ks exposure taken between Dec. 
12, 1992 - Jun. 14, 1993; Wang 1995; Wang & 
Gotthelf 1998), rh400779n00 (26 ks between Aug. 
21, 1996 - Jan. 27, 1997) and rh400779a01 (79 ks 
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between April 3-29, 1997). The first two expo- 
sures are, however, too short to be effective for 
detecting the pulsar. Our timing analysis was 
thus based exclusively on the latest observation. 
During the observation, the satellite was wob- 
bled, and PSR J0537-6910 was located very close 
(~ 2'.7) to the time-averaged optical axis. The 
expected Point Spread Function (PSF) is about 
6" (FWHM). 

We processed the data to optimize their tim- 
ing and positioning accuracy. First, we deleted 
time intervals with unstable detector roll angles 
and/or with fewer than 3 guide stars for posi- 
tion tracking, and shifted counts in individual 
observing segments (typically about 1 ks long) 
to have a common centroid position of the com- 
pact source (Fig. 1). This process improved 
the photon position accuracy by about l" in 
terms of the FWHM of the PSF. The source 
is not strong enough, however, to allow for a 
wobble phase-resolved centroid correction; thus 
significant wobble phase-dependent pointing er- 
rors may still be present in the data. Second, 
we corrected the absolute astrometry of the data 
by matching the X-ray and optical positions of 
three identified objects (the foreground stars cal 
69 and cal 71 as well as the LMC 50 ms pulsar 
B0540-69; Wang 1995). After this correction, the 
residual differences between the X-ray and opti- 
cal positions of these three objects are all less 
than ~ C//4. Third, to reduce the non-cosmic 
contamination we selected data only within the 
pulse height channel 2-10. Fourth, we converted 
the arrival time of each photon to the solar sys- 
tem barycenter using the JPL DE200 ephemeris. 

3. Data Analysis and Results 

For the timing analysis, we extracted counts 
within a circle of 7" radius around the centroid 
of the compact source — the suspected location 
of the pulsar (Wang & Gotthelf 1998). These 
counts were to be used to search for the periodic 
signal around the pulse frequency (y) extrapo- 
lated from the previously measured values and 



the mean frequency derivative i> = —1.9746 x 
10" 10 s" 2 (Marshall et al. 1998). However, we 
found that a small uncertainty in z> could cause a 
significant phase shift over the observation's to- 
tal duration of 2234676 s, during which individ- 
ual observing segments were separated by vari- 
ous data gaps resulting from the satellite passage 
through the South Atlantic Anomaly, Earth oc- 
cultation, etc. Based on the scatter of measured 
pulse frequencies around the linear spindown re- 
lation (Marshall et al. 1998), we estimated that 
the v uncertainty is ~ 0.04%, which could cause 
a phase shift of 0.2 cycle. Thus we removed sev- 
eral small observing segments (total 266 counts) 
that were largely separated from the main mid- 
dle portion of the observation. The remaining 
1288 counts are all within MJD 50545.0-50554.8 
(849270 s). Over this time interval the phase 
shift of 0.03 cycle is sufficiently small, compared 
to the pulse width of ~ 0.1 cycle (FWHM). 

We then constructed a frequency-gram (x 2 vs. 
u) of the selected counts to search for the peri- 
odic signal. For each trial frequency, we folded 
the counts into 10 phase bins to optimize the 
signal in the counting noise-dominated data and 
calculated the x 2 statistic of the resultant pro- 
file to test against the null hypothesis of no sig- 
nal. A significant signal with x 2 ~ 60 (^ 12cr) 
is clearly present in the frequency-gram. Due 
to the presence of noise in the data, the peak 
with the maximum x 2 value, however, is not nec- 
essarily the optimal measurement of the signal 
period. In fact, the exact peak period depends 
slightly on the frequency step used to construct 
the frequency-gram. Both the noise and the fre- 
quency step dependency can be suppressed by 
convolving the frequency-gram with an appropri- 
ate filter. We generated such a filter directly with 
the smoothed signal peak in the frequency-gram 
with the smallest frequency step. Because the 
shape of the peak depends primarily on the pulse 
profile and on the total observing duration, the 
filter can be re-binned for frequency-grams with 
larger steps. We adopted the signal peak centroid 
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of the filtered frequency-gram as our period mea- 
surement (Table 1), which is nearly independent 
of its chosen frequency step and is insensitive to 
the exact shape of the filter. An analog of this 
approach is the position centroid determination 
of a source in a 2-D X-ray image convolved with 
its PSF, or simply with a Gaussian filter. 

The uncertainty in our period measurement 
of the pulsar is determined through Monte Carlo 
simulations, based on the measured pulse pro- 
file before the binning and on the actual observ- 
ing segments. Each of our 1000 simulated data 
sets had the same number of counts as in the 
observed data, and were analyzed in the exactly 
same way. The dispersion in the obtained pe- 
riod distribution was adopted as an estimate of 
the uncertainty in our measurement of the pulsar 
period (Table 1). 

Table 1: Pulsar Parameters 
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Fig. 2. — Pulse profile of PSR J0537-6910 in the 0.1- 
2 keV band from the RHRI observation. The dashed line 
refers to the mean number of counts per bin averaged 
outside the 0.55-0.8 on-pulse phase interval. For a bet- 
ter illustration of the profile in this plot, a full cycle is 
duplicated and 25 bins per cycle are used. 



R.A. (J2000) 
Dec. (J2000) 

Frequency v (MJD 50540.5) 
Mean v (1993-1997) 
or Mean P 

Pulse Width (FWHM) 
RHRI Count Rate 
Pulsed Flux (0.1-2 keV) a 



5 A 37 m 47 s 2 

-69°10'23" 
62.05213169(±8) s" 1 
-1.97493 x 10~ 10 s~ 2 
5.1271 x 10~ 14 s s" 1 
~ 2 ms 

3.0 x 10~ 3 counts s" 1 
~6x 10~ 13 
ergs cm -2 s _1 



"Assuming a power law spectrum of energy slope a = 0.6 
and an absorption of 1 x 10 22 cm -2 (40% solar abundance; 
Marshall et al. 1998.; Wang & Gotthelf 1998). 

Fig. 2 shows the pulse profile. The pulse pro- 
file is narrow (FWHM ~ 0.1 duty cycle) and re- 
sembles the ones seen in the previous detections 
(§1), all of which are in energy bands ;> 2 keV. 
The pulsed intensity is about 30% of the total in- 
tensity during the full-width duty cycle of ~ 0.25, 
or about 10% if averaged over the full cycle. Fig. 
2 also indicates the presence of a possible weak 



interpulse, about 0.5 phase-shifted from the main 
pulse; but the significance of the detection is low 
(~ 3a). 

There is little doubt that the pulsar is in- 
deed within our adopted count collecting aper- 
ture. We find that the pulsed signal diminishes 
with the increase of the aperture radius, appar- 
ently due to the proportional increase of un- 
pulsed emission from the remnant. To determine 
the pulsar position accurately, we produced an 
image corresponding to the pulsed emission only. 
The insert of Fig. 3 includes an on-pulse intensity 
image with the off-pulse data subtracted. The in- 
tensity was calculated within the phase interval 
between 0.55 and 0.80. The average radial in- 
tensity distribution around the source centroid 
in the image is consistent with the PSF (Fig. 3). 
We thus adopt the source centroid as the posi- 
tion of the pulsar. The statistical uncertainty in 
the position is ~ 1". But we could not rule out a 
systematic error up to ~ 2", introduced by such 
effects as the asymmetric shape of off-axis PSF, 
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which may affect the accuracy of our astrometry 
correction (§2). Fig. 3 also includes an off-pulsed 
intensity image which clearly shows a strong in- 
tensity peak, probably barely resolved, and an 
elongated feature that extends at least ~ 30" to 
the northwest. 

4. Discussions and Conclusions 

The above results, as summarized in Table 1, 
provide an accurate position of PSR J0537-6910 
inside the SNR N157B, represent the first detec- 
tion of the pulsar in the 0.1-2 keV band, and 
extend the period measurement of the pulsar 
to Epoch 1997. This latest measurement, to- 
gether with the earliest ASCA observation (June 
13, 1993; Marshall et al. 1998), gives an aver- 
age period derivative (Table 1) that is slightly 
greater than 5.126 x 10~ 14 s s _1 as derived from 
the ASCA and 1996 RXTE observations (Mar- 
shall et al. 1998). This discrepancy is consis- 
tent with the presence of a giant pulsar glitch 
(8v/v ~ 10~ 6 ) occurred between the latter pair 
of observations. The pulsed emission accounts for 
only ~ 3% of the total luminosity of SNR N157B 
in the RHRI 0.1-2 keV band (Wang & Gotthelf 
1998). This pulsed contribution is substantially 
smaller than ~ 10% in the ASCA (2 - 10 keV) 
band, but is consistent with the extrapolated 0.1- 
2 keV flux from the high energy band (Table 1). 
The pulsed emission spectrum (power law energy 
slope a « 0.6) is significantly harder than the 
mean remnant spectrum (a ~ 1.5; Wang & Got- 
thelf 1998). The 0.1-10 keV luminosity of the 
pulsed emission (into Air) is ~ 3.3 x 10 35 ergs s _1 , 
or ~ 4 x 10~ 4 of the pulsar's spin-down luminos- 
ity. 

These new results are fully consistent with 
our early suggestion about the dynamic relation- 
ship between the pulsar and the Crab-like rem- 
nant N157B (Wang & Gotthelf 1998). First, 
as indicated by a separation of ~ 12" between 
the X-ray and radio emission peaks of the rem- 
nant, the pulsar is likely moving at a high ve- 
locity (~ 10 3 km s _1 ). Second, this motion has 



led to the formation of a bow shock around 
the pulsar. The enhanced synchrotron radia- 
tion of relativistic particles in the shocked pul- 
sar wind accounts for the unpulscd component 
of the compact source. Third, a collimated out- 
flow, driven by the pressure difference between 
the bow shock and the ambient medium (Wang 
& Gotthelf 1998; see also Wang et al. 1993), pow- 
ers an extended pulsar wind nebula where most 
of the pulsar spindown energy has been dumped. 
The separation between the pulsar and the main 
pulsar wind nebula explains the elongated fea- 
ture seen in both radio and X-ray. While both 
the pulsar position and the small pulsed con- 
tribution to the compact source are confirmed 
here, we continue to study the relationship be- 
tween PSR J0537-6910 and N157B using archival 
data in various energy bands and with upcoming 
AXAF observations. 
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Fig. 3. — The pulsed and unpulsed soft X-ray emissions 
from N157B. The pulsed emission was derived from the 
RHRI data within the full-width of 0.25 duty cycle and 
after subtracting normalized off-pulsed data. The average 
pulsed and unpulsed radial intensity profiles are marked 
with crosses and triangles, respectively. The pulsed emis- 
sion data are reasonably well-fitted (x 2 /d.o.f. = 24T./22) 
with the PSF shown as the solid curve. The slight excess 
of the data points above the PSF between 5-10" can be 
explained by telescope wobble-related aspect errors (§2). 
For comparison, the dotted curve shows the PSF normal- 
ized to the intensity of the unpulsed intensity profile. The 
inserted plot shows the intensity images of pulsed (a) and 
unpulsed (b) emissions. The contours are at 3, 6, 12, 24, 
and 48 xlO -2 counts s _1 arcmin~ 2 . 
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